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反对 称 铺设 复合 材料 层 合板 非 线 性 后 屈曲 分 析 
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Ja 要 :基于 层 合 板 党 理论 ,考虑 反对 称 铺设 层 合板 的 拉 谊 耦合 效应 和 后 必 曲 过 程 中 的 非 线 性 几何 
变形 ,推导 了 由 应 力 函 数 和 挠 度 表示 的 复合 材料 层 合 板 的 后 属 曲 控制 方程 。 引 入 无 量 纲 参 数 对 控 
制 方程 和 边界 条 件 进 行 无 量 纲 化 ,以 消除 材料 参数 及 几何 尺寸 对 分 析 结 果 的 影响 。 采 用 摄 动 法 将 
无 量 纲 的 非 线 性 控制 方程 及 边界 条 件 展 开 成 一 系列 非 齐 次 线性 摄 动 方程 组 ,分 析 各 阶 摄 动 方程 的 
通 解 与 特 解 的 构造 ,并 逐次 求解 ,建立 了 反对 称 铺设 复合 材料 层 合板 受 单 向 均 布 压力 作用 的 临界 层 
曲 荷载 及 后 慑 曲 平 衡 路 径 的 理论 解 。 进 而 运用 ABAQUS 软件 对 复合 材料 层 合板 在 面 内 压缩 载荷 作 
用 下 的 属 曲 和 后 属 曲 进行 有 限 元 分 析 , 结 果 表 明理 论 解 与 ABAQUS 结果 十 分 接近 ,验证 了 理论 解 的 
正确 性 。 在 此 基础 上 进一步 讨论 了 铺设 角度 、 铺 设 层 数 和 拉 谊 耦合 效应 等 对 层 合 板 后 属 曲 性 能 的 
影响 。 研 究 发 现 层 合板 的 属 曲 载荷 受 铺设 角度 与 层 数 的 影响 较为 显著 ,而 拉 谊 耦合 效应 使 板 的 必 


后 强度 大 大 降低 。 
关键 词 :复合 材料 层 合板 ;后 慑 曲 ; 摄 动 法 ;有 限 元 分 析 
中 图 分 类 号 :0343.5 ”文献 标志 码 :A DOI:10. 11776/j. issn. 1000-4939. 2023. 04. 006 


Nonlinear post-buckling analysis of composite laminates 
with antisymmetric angle-ply 
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(MOE Key Laboratory of Disaster Forecast and Control in Engineering ,School of Mechanics and Construction 


Engineering, Jinan University ,510632 Guangzhou , China) 


Abstract: Based on the theory for laminated plates and shells , governing equations for post-buckling of an- 
ti-symmetrical laminated composite plates are derived in terms of stress function and transverse deflection 
by considering the tension-bending coupling effect and the von-Kaman nonlinear geometric relation. By in- 
troducing the non-dimensional parameters ,the nonlinear governing equations and the boundary conditions 
are nondimensionalized to eliminate the influence of the material properties and the geometric parameters. 
Perturbation approach is applied to expand the governing equations as well as the boundary conditions into 
a series equations of a small perturbation parameter. Analytical solutions of the critical buckling load and 


the post-buckling equilibrium path for anti-symmetric cross plied composite laminates with simply supported 
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boundary conditions are obtained by examining the characteristics of the general and particular solutions. 


After that ,the buckling and post-buckling of the composite laminates with antisymmetric angle-ply under 


in-plane compressive load are analyzed using ABAQUS software, and the results from ABAQUS are 


compared with the theoretical solutions. It shows that the two results are in very good agreement which 


verifies the accuracy and effectiveness of the proposed analytical approach. Using the established 


solutions ,the influences of the orientation angle of the composite ply ,the number of the composite ply and 


tension-bending coupling effect on the post-buckling performance of the laminate are discussed. It is found 


that the critical buckling load of the laminate plate is significantly affected by the orientation angle and the 


tension-bending coupling effect will greatly reduce the post-buckling load carry capability of the laminates. 


Key words :composite laminate; post-buckling; perturbation approach; finite element analysis 


复合 材料 层 合 板 壳 结 构 ,因为 轻 量化 .高 比 强度 
和 比 模 量 等 优点 ,通常 作为 承载 部 件 广泛 的 应 用 于 
机 械 化 工 .航空 航天 等 领域 。 承 受 面 内 压缩 荷载 作 
用 的 板 壳 结 构 , 当 承载 力 增加 到 某 一 极限 载荷 值 时 ， 
结构 很 容易 会 发 生 失 稳 破 坏 ,进而 给 工业 和 日 常 使 
用 带 来 极 大 危害 , 因此 层 合 结构 的 稳定 性 分 析 极 具 
工程 价值 。 

20 世纪 初 Von Kármán 建立 了 著名 的 板 壳 非 线 
性 弯曲 理论 '" ,奠定 了 板 壳 几何 非 线 性 分 析 的 基础 ， 
从 此 世界 各 国学 者 们 对 复合 材料 层 合板 壳 的 非 线性 
力学 性 能 开始 进行 深入 研究 Б ЈЕ 
壳 大 拨 度 弯曲 时 引起 的 中 面 应 变 , 视 双 模 量 扁 过 为 
两 种 材料 组 成 的 层 合 扁 壳 ,解决 了 双 模 量 简 支 扁 这 
在 均匀 内 压 作 用 下 的 非 线性 弯曲 问题 。BELLIFA 
等 5 根据 Eringen 的 非 局 部 弹性 理论 ,考虑 剪 切 力 的 
轴 向 位 移 效 应 ,提出 了 一 种 针对 纳米 尺度 梁 的 非 线 
性 后 届 曲 行为 的 非 局 域 零 阶 剪 切 变形 理论 。 张 作 亮 
等 利用 Von Kármán 大 挠 度 的 理论 ,结合 达 朗 贝尔 
原理 和 薄板 振动 理论 ,建立 了 爹 形 张 拉 膜 结构 在 冰 
和 置 丛 载 作用 下 的 非 线 性 动力 响应 控制 方程 。 当 板 这 
结构 处 于 小 变形 和 适度 旋转 的 情况 下 时 ,如 初始 后 
届 曲 阶段 ,采用 Von Kármán 非 线 性 几何 关系 可 以 获 
得 有 效 的 结果 ,但 当 板 壳 结 构 处 于 深 后 届 曲 阶段 时 ， 
则 必须 考虑 板 壳 结 构 的 大 变形 。SHARIYAT 和 
АЅЕМІ 利用 全 格林 应 变 张 量 来 描述 板 壳 结构 深 后 
届 曲 过 程 中 的 几何 非 线 性 关系 ,并 分 析 了 受 单 轴 和 
双 轴 压缩 载荷 作用 下 变 刚 度 矩 形 板 的 后 屈曲 行为 。 


Mindlin 模型 ,同样 利用 有 限 元 法 分 析 了 含有 圆 形 和 
椭圆 形 分 层 的 复合 材料 层 合 板 后 屈曲 行为 。 卢 讲 瑞 
等 "用 有 限 元 法 研究 了 关于 无 人 机 3 个 关键 部 位 
的 复合 材料 蒙 皮 壁 加 筋 板 的 非 线 性 届 曲 问题 ,得 到 
了 壁 板 的 屈曲 特性 和 后 屈曲 损伤 演变 过 程 。 林 国 伟 
等 "和 和 孙 中 雷 等 '” 均 采用 了 渐进 损伤 分 析 方 法 ， 
分 别 对 复合 材料 加 筋 板 在 承受 轴 压 载荷 下 的 届 曲 载 
荷 及 其 破坏 模式 和 后 届 曲 极限 破坏 载荷 及 其 失效 模 
式 进行 了 深入 的 研究 ,前 者 改进 了 已 有 的 工程 计算 
方法 ,使 计算 结果 与 实验 值 更 加 吻合 "i ,后 者 分 析 
了 各 种 参数 对 复合 材料 加 筋 板 后 届 曲 特性 的 影 
HU ,为 复合 材料 加 筋 板 的 设计 提供 了 参考 依据 。 

在 分 析 结 构 承 受 复杂 外 力作 用 下 的 力学 性 能 
时 ,有 限 元 法 具有 不 可 替代 的 优势 ,而 近似 法 ,如 作 
辽 金 法 ' ”| Koiter 渐进 法 '"“ 、 摄 动 法 等 , 则 能 更 深 
入 解释 结构 失效 的 机 制 ,揭示 结构 内 部 的 受 力 情况 。 
摄 动 法 是 当前 解决 工程 技术 问题 和 科学 问题 的 主要 
数学 工具 之 一 ,因此 被 广泛 的 用 于 求解 结构 的 非 线 
性 问题 。 祭 桂林 '" 基于 Euler-Bernoulli 梁 理论 研究 
了 热膨胀 下 梁 的 横向 非 线 性 振动 问题 ,并 采用 二 次 
摄 动 法 进行 求解 ,与 传统 摄 动 法 的 计算 结果 进行 了 
对 比 验证 。 李 萍 等 "采用 小 参数 援 动 法 研究 了 脱 
层 对 复合 材料 层 合 梁 的 后 届 曲 的 影响 。YANGA J 
等 ("采用 摄 动 法 讨论 了 功能 梯度 和 矩形 板 的 大 变形 
问题 ,再 基于 四 边 固 支 的 边界 条 件 ,讨论 了 在 预 加 横 
向 载荷 和 面 内 压缩 载荷 作用 下 的 后 屈曲 行为 。 
DASH 等 ”| 基于 拉 格 朗 日 非 线 性 应 变 -位移 关 系 的 


由 于 结构 非 线 性 弯曲 问题 的 控制 方程 难以 获得 
精确 解 ,学 术 界 常常 采用 数值 法 ,如 有 限 元 法 、 渐 i 
损伤 分 析 法 ,来 进行 分 析 。ALMITANI' 采用 有 限 
分 析 方法 研究 了 对 称 与 反对 称 功能 梯度 梁 的 屈曲 | 
题 。 杨刚 等 1 提出 了 改进 层 合板 壳 的 二 维 分 层 


ЫЙ: 


2] 


高 阶 剪 切 变形 理论 ,利用 摄 动 法 研究 了 层 压 复合 板 
的 屈曲 和 后 屈曲 问题 。 夏 飞 等 ”同样 运用 摄 动 法 
研究 了 在 横向 剪 切 与 湿热 环境 下 温度 、 脱 层 条 件 等 
各 参数 对 复合 材料 层 合板 的 届 曲 与 后 届 曲 的 影响 。 

由 于 正 交 对 称 铺设 的 层 合板 不 存在 拉 弯 耦合 效 


投稿 网 站 :http://ejam. хўш. ейи. en 微 信 公众 号 :应 用 力学 学 报 


780 应 用 力学 学 


应 ,因此 ,上 述 文献 在 使 用 摄 动 法 时 得 到 的 非 齐 次 线 
性 摄 动 方 程 组 ,可 以 通过 对 应 力 函 数 与 找 度 进行 解 
耦 来 求解 。 而 当铺 设 情 况 为 反对 称 铺设 时 , 拉 弯 耘 
合 效应 的 存在 导致 应 力 函 数 与 挠 度 无 法 解 耦 。 本 研 
究 在 轴 压 作用 下 的 简 支 反对 称 铺设 复合 材料 矩形 层 
合板 的 后 届 曲 问题 ,将 Von Kármán 非 线性 几何 方程 
引入 到 层 合板 的 平衡 微分 方程 中 ,建立 了 反对 称 铺 
设 层 合板 的 非 线 性 后 届 曲 控制 方程 。 通 过 使 用 二 阶 
小 参数 摄 动 法 ,求解 非 线性 控制 方程 组 ,得 到 了 临界 
届 曲 载 集 与 后 届 曲 平衡 路 径 。 应 用 ABAQUS 软件 
进行 仿真 分 析 , 并 将 有 限 元 解 与 理论 解 进行 比较 , 验 
证 理论 解 的 准确 性 与 可 靠 性 。 最 后 进行 算 例 分 析 ， 
讨论 不 同 铺设 方式 、 拉 弯 耦 合 效应 等 对 层 合板 的 临 
届 曲 载荷 及 后 届 曲 性 能 的 影响 。 


1 基本 方程 


1.1 模型 的 建立 


考虑 图 1 所 示 的 四 边 简 支 的 复合 材料 层 合板 ， 
长 度 为 a, 宽度 为 5, 铺 层 数 为 N, 总 厚 为 i, 单 层 板 厚 
为 h, 受 x 方向 均 布 荷 载 NN, 作用 ,建立 直角 坐标 系 。 
各 子 板 的 几何 尺寸 如 图 1 所 示 。 


(b) 层 合板 的 载荷 位 移 边 界 示意 图 
图 1 层 合板 模型 示意 图 


Fig.1 Model of rectangular plate 


12 层 合 板 的 非 线性 大 挠 度 方程 


设 в, лет ew 为 层 合板 的 中 面 应 变 , к, K, .Kw 为 
中 面 的 曲率 、 扭 曲率 ,uwv、w а i 
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则 由 Von Kármán 非 线性 板 壳 理论 可 得 


в, u, + бр 0?/2 
- | І p Ë Са 22/2 | - 
e, u, + Ç, + W W, 


K, =w ‚хх 
H: Е си a 
к -w , 


式 中 ,uw KIEM и 对 变量 x Ае. нане, 从 
式 (1) 中 消去 4 和 vw, 得 到 层 合板 非 线 性 弯曲 的 协调 
方程 为 


£, Të, 2 Eyo = w? 7 


- 00у (3) 
本 研究 反对 称 铺设 各 向 异 1 性 层 合板 , 设 0 为 第 
层 单 层 板 的 坐标 系 与 材料 主 方向 的 夹 角 , 则 该 层 的 


弹性 刚度 为 


Го, 
A r ОС s 465? 
le съз S лее oh 
Е 
Q> Ss 2с? С 4С°8° 0, 
0 | cas се -Cs -CS —2С8(С° - 5°) | Q, 
16 
Е CS (8-С -CS 2CS(C -S) LO, 
О» Wasa -D0282 CS (C – 5)? 
L Q% 
(4) 
式 中 :C = cos0 ;S = sin0 , H. 
2 Е, = E, _ с; 
É = 8 = Ts = 12 
Шоо] Шоо 
О = Q = Кл Qn ‚Мо = Wabi 
(5) 


式 中 :Ei、E, 表 示 材 料 在 1.2 弹性 主 方向 上 的 弹性 模 
量 ;Co 表示 材料 在 12 平面 内 的 剪 切 弹性 模 量 ;jw 
和 pw 为 不 同方 向 的 泊 松 比 。 

根据 复合 材料 力学 理论 2 
程 为 


, 层 合板 的 本 构 方 


A. b. 0 
ls ol ® 
式 中 : 为 层 合 板 的 中 面 薄 膜 内 力 ; 为 层 合板 的 内 
力矩 ;4; 为 拉 伸 刚度 ; B, 为 耦合 刚度 ; D; 为 弯曲 刚 
度 ,分 别 表示 为 
力学 学 报 


ui 
үт 

心 
ж 


> Qp 2 – 


ï T у (Q), - 2.) (7) 


оо 
1 


Р; = т ry (QD 12 a) 


式 中 ,za 表示 第 层 单 层 板 。 根 据 经 典 板 党 理论 , 层 
合板 后 屈曲 的 平衡 方程 为 


Nu +N,,,=0 
N, +N,,=0 
M... +2M ¿+M „+ „+ № „+2ХМ мю =0 
(8) 
引入 应 力 函 数 
N =P Nb u Ng =D „. (9) 


Уу? 


将 式 (6) 和 (9) 代入 到 协调 方程 (3) 和 平衡 方程 
(8) 的 第 3 式 中 ,同时 考虑 到 对 反对 称 铺设 层 合板 ， 
有 Bi =b, = В = В, =0 
得 到 反对 称 铺设 层 合 板 非 线性 后 届 曲 方程 为 


D u 


АР, Рю у + (2р +4020 „уу 一 


(ОВ, -BOB - (ОВ, - Bp) P ay- 


P „w - b awp +2Ф p way =0 
(10) 
А» DB + Ar Ф „„ + (2А + А „)Ф ‚ххуу 
(2 Bie В) хууу + (2 B, – В, )%0 ххху 
2 
W ху š W ,er yy 
(11) 
式 中 ,有 =B4- ,В-А'В,р = -BA B +D, 
本 研究 对 象 为 四 边 简 支 的 层 合 板 , 边 界 条 件 可 
表示 为 
w 三 
* =0, a:f I (12) 
M. =В„Ф „+ „ +D u „у =0 
№ = 
y=0, | | | (13) 
М, =В»Ф y +0 № a + Dw „у =0 


2 摄 动 法 求解 


本 研究 采用 摄 动 法 对 式 (10) 和 (11) 进行 求解 。 
引入 无 量 纲 参数 , 即 
=t/a,B=t/b,w = w/t, 


x=x/a,y=y/b,a 
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D; = Р„/Аь? A. =А„!/А„ .B =B; = BZ, 
В, = B ,/t,$ = Ф/Аь? ,N, = N.P /А„? (14) 
НЯ Оо, Ф, М, 展开 式 为 
W =W W, Wg + шуш шр + изии + `° + WW" Ww, 
Ф=Ф, + + Фүш ш, + B ww + Dyw ws + 

+ Фб, 
N. =N. +N nw w +N шш + ` + Муш” ш?" 


(15) 

式 中 :wo = (w) Zt 为 摄 动 小 参数 ; N。 是 反对 称 铺设 
复合 材料 层 合板 量 纲 为 1 的 屈曲 荷载 。 

将 上 述 无 量 纲 函数 展开 式 代 入 控制 方程 ,再 对 
小 参数 wo 的 不 同 次 窜 进 行 合 并 ,得 到 各 阶 摄 动 方 程 
组 ,代入 边界 条 件 , 逐 阶 求解 便 可 获得 足够 精确 的 
解 。 附 录 A 给 出 了 高 阶 求解 结果 。 

零 阶 方程 OG) 为 
(2 Bo - Вь)о?В $... +(2 В, -В,,)ов? Ф 
А» о“ Ф, „а +А 8° Ф + (2 Ар + 
Аб, PB Фу ы, =0 


0,хууу =0 , 


0,yyyy 


(16) 
满足 边界 均 压 应 力 为 
P, = -М„у?/2. (17) 
一 阶 方程 0(w,o ) 为 
Draw + Dab w, „уу, + (20, +400) PBW y = 
(2. Bo -Вь)о?В Ф, - (2Вь - Ba) aB Ф, – 
cg yy w, ax T Ф, И? РЭ +208? Ф, Елее =0, 
Аза Bi wn +А В" Ф, »» + CAR + Е 
Ав PP Bi + (2 В, — Вь) ав ww + (2 B, - 
Ba) BW, „„=0 


(18) 
满足 边界 条 件 的 解 为 
wi =w,sin( mT x) sin(nm y) 
Ф, = ф,ю„соз( тт х) соѕ(пт у) (19) 


其 中 
(2 B, А: Boe) aB mn + (2 В» = В, а Bmn 


Pa = — = =E 
Azam? + A B ni (2А + Ad )aB mn 
(20) 
Б p 2 2 р. 2 4 _ _ 

N. = т a + ü + (2D +4р,„)п| + 

3 

атп, П, 

T| OBa - B.) B T (2B。 Е В») Ве, 


(21) 
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二 阶 方程 0(wi) 为 
Бо ш, „ы, + Ре + 2D, +4, За B Ww Е 


(2B。 F B,D BD, „, Е OBa Е BB D, „, Е 
apB Ф, w, a + Ф, Wi a + Ф, 


2 ВФ 
А» а*Ф, o + АВФ, „„ + OQA} + ADEE D, „ + 


w, „+ 
уу 


0,уу 1,yy 1,хх 


ш „ = 0, 


І, ху 


(2 В; Е Bo ) opB Ww, „, + (2 B, Бы В, 20 Во, s = 
opB Cw?, = ПОЛ 


(22) 


满足 边界 条 件 的 解 为 
w, = 0, 
Ф, = paw соѕ(2тлтх) + рю сов(2тт у) - 
№? 7 /2 (23) 

其 中 I 

Фи = B?n'/32A;1 om, 

Фф» = оёт? /З2А Вп? (24) 

三 阶 方程 0(wi) 为 

Юс юу ы + РВ yy + (2р, +4D,)o Ваз, = 
(2Bo = B,)a°8 Фу м, T (2В,, = By ) ap’ Ф, on Е 
ов? ( Ф, W3 a + $, , w, T Ф, „w, ж) =0, 
Az a By FAD B Ds „+ (245 + 
Ав, ) PP D, „„ +(2 В, Е Bo) aB ws oyy (2 В» С 
Bu ) o° Du, ы, =0 


3 ,xxxx 


(25) 
满足 边界 条 件 的 解 为 
w, =шуүюувїп(3ттх)зїп(пту) + 
wyw sin( тлтх) sin( 3nqry) ， 
Ф, = фути? соѕ(Зтлтх) cos(nmy) + 
gaw соѕ( тях) cos (Злу) (26) 


同 理 ,四 阶 摄 动 方程 0(wi ) 为 
Da Q W4 м + DaB w „уу + (2D +40,)о Вад у = 
(2В,, - Bio) a BD; ay - (2В,, - Ba) AB D „‚- 
ap ( Do ws, +Ф 
208 ( Ф 
Az oa By „ ЖАЛ В'Ф, + (245 +AG )o ВФ, t 
(2 Bis — Bo) op w a, + (2 В, -Ba )a B w, = 


252 Fa ысы оа 
а В ( =W] „103 уу — W3 Wi yy) 


1,yy W3 ax + Ф, шу + Ф, „03у F 


W3 xy + P; Wi ay) =0, 


1›ху 


(27) 
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满足 边界 条 件 的 解 为 

w, = wawisin(4mmx)sin(2nmy) + 
wpw sin(2mnx)sin(4nmy) + 
10450 5іп(2тлтх) sin(2nmy) 


Ф, = фии соѕ(Аттх) сов(2пату) + 


Фош соѕ(2тлтх) соѕ(Апту) + 

Фази соѕ(2тлтх) соѕ(2пту) + 

Paw соѕ(Атлатх) + фахш соѕв(Апату) + 

фаи соѕ(2тлтх) + фази соѕ(2лату) -Naw у 2/2 
(28) 


至 此 ,已 经 可 以 写 出 大 挠 度 渐 近 解 
ю = wolw,sin(mmx)sin(nmy)]1+ 


5 ww віп (Зтттх) sin( пту) + 


0 2 = 
wyw sin( тттх) sin( 3nqry) 


10410 ( вїп4#тттх) вїп(2пту) + 
w wow sin(2mmx ) sin(4nmy) + | 十 Ось) ， 


ww віп (2тлтх) sin(2nmy) 


S| 
I 


— N.” /2 + ше! ло, cos( тах) cos( ny) 1+ 
1 w соз(2ттх) + pyw cos(2nmy) ] 

№?у/2 

3 


a [ paw cos(3mmx)cos(nmy) + 
“| |» 


gaw соѕ( тах) сов(3птту) 

wil фа шсоѕ( Атах) со5(2пту) + 
Фош соѕ(2тлх) соѕ(Апту) + 

Фази соѕ(2тлтх) соѕ(2пту) + 

Фаш соѕ(Атлтх) + paswicos(4nmy) + 


Флеш соѕ(2тлтх) + gywtcos( 2nay) 一 


Мышу /2] + Ос) (29) 
ER 0 (х,у) = (1/2т,1/2п) ,有 
Wn = W,Wo 一 Сз + Wwa VW W (30) 


IN (30) 中 w, w 可 视 为 w, 的 表达 式 , 即 
WaWo = Wp + Сю, + wW, 
代入 式 (15) 第 3 式 ,可 得 到 板 的 后 屈曲 的 平衡 路 径 
N. =N. +N +[N+2N( +w) Jw, + 
(31) 


з 计算 与 结果 分 析 


本 研究 选取 四 边 简 支 的 Т300/5208 石墨 环 氧 复 
合 材 料 层 合板 作为 研究 对 象 , 层 合板 的 几何 参数 为 : 
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a=b=2 m,h =0. 000 138 m; 材 料 参数 为 :已 ，= 
9 GPa, E, =360 GPa, Gp =4.5 CPa,C =1.8 GPa, 
Gis = Ga, U =0.25,и„ = и,Е,,/Е,, о 

为 了 验证 分 析 的 准确 性 ,运用 ABAQUS 有 限 元 
软件 建立 模型 并 进行 分 析 。 在 屈曲 分 析 步 中 , 边界 
条 件 设置 为 四 边 简 支 ,在 层 合板 模型 x 轴 向 的 两 边 ， 
设置 1 N/m 的 压缩 载荷 ,最 后 输出 的 特征 值 即 为 层 
合板 的 届 曲 载 答 值 ,线性 摄 动 的 最 大 迭代 次 数 设 为 
1 000 次 ,以 确保 模型 计算 结果 的 收敛 性 。 

层 合板 的 后 屈曲 仿真 分 析 是 建立 在 已 求 得 屈曲 
载荷 和 届 曲 模 态 的 基础 上 进行 的 ,采用 静态 Riks 法 
分 析 层 合板 的 非 线性 后 屈曲 问题 。 与 线性 摄 动 一 届 
曲 分 析 不 同 ,在 静态 Riks 分 析 步 中 , 层 合板 必须 要 
有 具 有 一 定 的 初始 缺陷 ,这 样 才 能 保证 层 合板 在 轴 疝 
压缩 荷载 的 作用 下 可 以 出 现 后 屈曲 挠 度 即 面 外 变 
形 。 在 已 有 的 有 限 元 屈曲 分 析 的 结果 基础 上 ,将 届 
曲 模 态 的 一 定 倍数 作为 初始 的 缺陷 赋予 到 后 屈曲 有 
限 元 计算 模型 ,然后 进行 加 载 分 析 。 同 时 ,在 Riks 
分 析 步 进行 计算 时 ,必须 考虑 几何 非 线性 选项 , 且 需 
要 设置 最 大 载荷 比例 系数 作为 终止 条 件 。 在 进行 计 
算 中 ,总 迭代 数 设 置 为 1 000 次 ,以 确保 获得 较为 完 
整 和 明显 的 后 屈曲 路 径 曲 线 。 在 层 合板 模型 x 方向 
两 边 施 加 线性 压缩 载荷 ,数值 设置 为 已 求 得 的 一 阶 
屈曲 载荷 值 ,进而 就 可 以 将 输出 的 板 的 挠 度 值 和 载 
荷 比 例 因 子 (LPF ) 进行 整理 ,绘制 在 不 同情 况 下 层 
合板 的 后 届 曲 路 径 曲 线 。 


3.1 有 限 元 验证 分 析 


图 2 给 出 了 N=6,0=45° 与 N=4,0=30° 时 后 
屈曲 阶段 载荷 - 挠 度 曲线 的 理论 解 与 有 限 元 解 的 比 
较 。 从 图 2 可 以 看 出 ,理论 分 析 与 有 限 元 分 析 预 测 
出 的 层 合 板 的 后 屈曲 路 径 是 一 致 的 ,两 者 都 表达 了 
以 轴 向 压缩 荷载 为 特征 的 后 屈曲 强化 特性 ,也 即 在 层 
合板 达到 屈曲 载荷 后 不 会 立刻 发 生 失 稳 , 而 是 随 着 面 
外 挠 度 变形 w 的 增加 ,外 载荷 也 在 以 某 个 幅度 增加 。 
两 种 方法 之 间 存 在 着 差 值 , 这 是 由 于 理论 计算 的 是 无 
初始 缺陷 的 平面 板 ,而 有 限 元 模拟 分 析 的 是 添加 了 初 
始 缺 陷 因 子 的 层 合 板 ,这 样 的 差 值 不 可 避免 。 

图 3 为 有 限 元 仿真 和 理论 分 析 计 算得 到 的 层 合 
板 后 屈曲 变形 模 态 图 ,图 中 算 例 选用 :V=4,0=45?"。 
图 3(a) 可 以 看 出 , 层 合板 在 达到 屈曲 后 ,将 发 生 面 外 
挠 度 变 形 ,呈现 的 是 整体 一 致 的 后 屈曲 变形 模 态 。 图 
3(b) 中 理论 计算 并 绘制 的 屈曲 变形 模 态 结果 与 图 3 
(a) 也 相同 ,再 次 验证 了 本 研究 理论 分 析 的 准确 性 。 
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(a) N=6, 0=45° 
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- - -MATLAB 
----ABAQUS 


(b) м4 0=30° 
2 有 限 元 结果 和 理论 解 得 到 的 后 屈曲 曲线 


Fie.2 Post-buckling curves obtained from finite 


element results and theoretical solutions 


U,Magnitude 
+1.000e+00 


+8.333e+02 
+0.000e+00 


(a) 有 限 元 仿真 解 
Е„/Е„=40 а=2. m,b=2 m,t=0.008 m 


00 
(b) 理论 解 
到 3 有限 元 结果 和 理论 解 得 到 的 屈曲 变形 模 态 


Fig.3 Post-buckline deformation modes of finite element 


results and theoretical solutions 
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3.2 算 例 分 析 


3.2.1 铺设 角度 对 属 曲 载荷 的 影响 
通过 有 限 元 软件 建立 了 复合 材料 层 合板 的 模型 
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合板 的 抗 弯 刚 度 增 大 ,临界 届 曲 载荷 呈 指 数 增长 。 
但 由 于 图 3 纵 坐标 为 无 量 纲 的 届 曲 载荷 ,这 使 得 当 
铺设 层 数 到 达 一 定 大 小 时 ,对 届 曲 载荷 的 增 大 不 再 
明显 ,可 以 看 出 有 限 元 分 析 结 果 与 理论 结果 吻合 ,两 


并 进行 有 限 元 分 析 来 验证 结果 的 准确 性 。 图 4 为 铺 
设 层 数 N =4,6,8 时 无 量 纲 届 曲 荷载 随 铺设 角度 变 
化 的 理论 解 与 有 限 元 结果 的 比较 。 可 以 看 出 ,当铺 
设 层 数 保持 不 变 时 ,理论 解 与 有 限 元 模拟 解 趋势 一 
致 。 当 铺设 角度 增 大 时 ,临界 届 曲 荷载 先 增 大 后 减 
小 ,在 90=4 "时 到 达 最 大 临界 届 曲 载荷 。 由 公式 
(21) 可知 ,临界 届 曲 载荷 主要 由 抗 弯 刚 度 决定 , 相 较 
于 其 他 铺设 角度 , 当 0 = 45 "时 , 层 合板 各 个 方向 的 搞 
SRBE D; 整 体 达到 最 大 ,使 得 临界 届 曲 载荷 最 大 。 


70 


60 


I= 50 


0 15 30 45 60 75 90 
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图 4 不 同 铺设 层 数 ,临界 屈曲 载荷 随 铺设 角度 的 变化 


Fig.4 Effect of fiber orientation angle on buckling load for 


different number of the laminate 
3.2.2 铺设 层 数 对 属 曲 载荷 的 影响 
图 5 显示 了 当铺 设 角 度 0 = 15°,45°,67.5° 时 ， 
无 量 纲 屈曲 荷载 随 铺设 层 数 的 变化 而 变化 。 
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图 5 不 同 铺设 角度 ,临界 屈曲 载荷 随 铺设 层 数 的 变化 


Fig.5 Effect of the number of the laminate on buckling load 


for different fiber orientation angle 


当铺 设 角 保持 度 恒 定时 , 随 着 铺设 层 数 增加 , 层 


者 误差 不 超过 3. 32% ,验证 了 理论 分 析 的 合理 性 。 
3.23 拉 谊 耦合 效应 对 后 属 曲 承载 能 力 的 影响 

为 了 研究 拉 弯 耦合 效应 的 影响 ,本 研究 分 析 了 当 
铺设 方式 分 别 为 [45 °/45 °] 对 称 铺设 和 [45 °/ —45 °] 
反对 称 铺设 时 ,正方 形 层 合板 的 后 届 曲 平衡 路 径 。 
图 6 给 出 了 两 种 铺设 情况 下 的 方形 板 的 后 届 曲 路 
径 。 图 6 结果 表明 , 拉 弯 耦合 效应 使 板 的 届 后 强度 
大 幅度 降低 。 
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a В=1.0, 0=45°, №2, В,=0 
ТЕЕ В=1.0, 0=45°, N=2 


0 (бэ Ш0 15 20 25 ЗО 35 0 
w 


6 PERA а Ji IH 3F-48F6 А W 
Fig.6 Influence of coupling effect of tensile on 
post-buckline equilibrium path 
3.24 铺设 角度 对 后 属 曲 承载 能 力 的 影响 
图 7 给 出 了 当铺 设 层 数 为 N=4 时 ,不 同 铺设 角 
度 的 复合 材料 层 合板 的 后 届 曲 平衡 路 径 ( 四 阶 近 
似 ) 。 


0 0.5 1.0 15 2.0 2.5 


87 铺设 角度 对 后 届 曲 平衡 路 径 的 影响 


Fig.7 Influence of fiber orientation angle on 


post-buckling equilibrium path 


由 图 7 А, 0 = 30° 0 =67. 5°48 ШК, 当铺 
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БЕШ 0 为 45° 时 , 层 合板 具有 较 高 的 临界 届 曲 载 答 ， 
但 后 届 曲 平衡 路 径 相 对 更 加 平缓 。 显 见 ,在 临界 届 
曲 载 荷 最 大 的 铺设 角度 时 反而 对 后 届 曲 的 承载 能 
相对 较 低 ,对 应 不 同 铺设 角 时 , 层 合板 后 届 曲 载 葆 一 
挠 度 曲线 具有 显著 差别 。 
3.2.5 铺设 层 数 对 后 屈曲 承载 能 力 的 影响 

图 8 显示 了 在 铺设 角度 96=45° 时 ,铺设 层 数 对 
后 屈曲 平衡 路 径 的 影响 。 可 以 看 出 当铺 层 数 V =2， 
4,6,8 时 , 铺 层 数 增加 ,导致 层 合板 刚度 增加 , 因此 
初始 届 曲 载荷 也 增加 ,而 随 着 铺 层 数 的 增加 , 拉 弯 
耦合 效应 减弱 ,因而 后 届 曲 平衡 路 径 逐 趋 平 缓 。 
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图 8 铺设 层 数 对 后 届 曲 平衡 路 径 的 影响 


Fig .8 Influence of the number of the laminate on 


post-buckling equilibrium path 
4 结 论 


本 研究 运用 复合 材料 层 合板 理 论 ,考虑 拉 弯 耦 
合 效 应 ,研究 了 反对 称 铺设 层 合板 的 非 线性 后 届 曲 
问题 ,通过 引入 Von Kármán 非 线性 几何 关系 ,建立 
了 层 合 板 的 非 线性 后 届 曲 控制 方程 ,并 应 用 二 次 摄 
动 法 进行 求解 ,获得 了 求解 反对 称 铺设 层 合板 的 临 
界 屈曲 荷载 及 后 届 曲 平衡 路 径 。 开 发 MATLAB $k 
件 进 行 参 数 分 析 ,详细 讨论 了 不 同 铺设 情况 bu gan 
合 效 应 等 对 层 合 板 的 届 曲 丛 载 与 后 届 曲 平衡 路 径 的 
影响 ,同时 将 理论 结果 与 有 限 元 结果 进行 对 比 ,验证 
了 本 研究 分 析 的 准确 性 与 合理 性 。 结 论 如 下 。 

1) 当铺 设 角度 增 大 时 ,临界 届 曲 丛 载 先 增 大 后 
减 小 ,在 0 =45 "时 到 达 最 大 临界 屈曲 载荷 。 铺 设 角 
的 变化 将 极 大 地 影响 板 的 届 后 强度 ,45° 为 最 佳 铺 
设 角 。 

2) 随 着 铺设 层 数 的 增 大 , 层 合板 的 临界 届 曲 载 
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荷 也 相应 地 增加 ,但 后 屈曲 强化 趋势 逐 趋 平缓 。 

3 ) 拉 弯 耦 合 效应 使 板 的 屈 后 强度 大 大 降低 , 增 
加 铺 层 数 可 减弱 拉 弯 耦合 效应 的 影响 。 

4) 与 9=30°* 和 0=67.5° 相 比 ,铺设 角 Ө H 45° 
时 , 层 合 板 具 有 较 高 的 临界 届 曲 载荷 ,但 后 届 曲 平衡 
路 径 相 对 更 加 平缓 , 届 后 强度 比 其 他 铺设 角度 的 强 
度 更 低 。 
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